library(tidyverse)
Registered S3 methods overwritten by 'dbplyr':
  method         from
  print.tbl_lazy     
  print.tbl_sql      
── Attaching packages ─────────────────────────────────────────────────────────────────────────────────────── tidyverse 1.3.1 ──
✓ ggplot2 3.3.5     ✓ purrr   0.3.4
✓ tibble  3.1.6     ✓ dplyr   1.0.6
✓ tidyr   1.1.3     ✓ stringr 1.4.0
✓ readr   1.4.0     ✓ forcats 0.5.1
── Conflicts ────────────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
x dplyr::filter() masks stats::filter()
x dplyr::lag()    masks stats::lag()
library(plotly)
Registered S3 method overwritten by 'data.table':
  method           from
  print.data.table     
Registered S3 methods overwritten by 'htmltools':
  method               from         
  print.html           tools:rstudio
  print.shiny.tag      tools:rstudio
  print.shiny.tag.list tools:rstudio
Registered S3 method overwritten by 'htmlwidgets':
  method           from         
  print.htmlwidget tools:rstudio

Attaching package: ‘plotly’

The following object is masked from ‘package:ggplot2’:

    last_plot

The following object is masked from ‘package:stats’:

    filter

The following object is masked from ‘package:graphics’:

    layout
library(rcartocolor)
merged <- read_tsv('../output/merged_table.tsv')

── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
  Specialty = col_character(),
  `No. of Programs` = col_double(),
  `Positions Offered` = col_number(),
  `Unfilled Programs` = col_double(),
  `MD Senior Applicants` = col_number(),
  `Total Applicants` = col_number(),
  `MD Senior Matches` = col_number(),
  `Total Matches` = col_number(),
  `MD Senior % Filled` = col_double(),
  `Total % Filled` = col_double(),
  `MD Senior Ranked Positions` = col_number(),
  `Total Ranked Positions` = col_character(),
  Year = col_double(),
  Class = col_character()
)
merged %>% 
  filter(Specialty == 'Dermatology') 
long_table_absolute <- merged %>% 
  select(Specialty,
         Year,
         Class,
         `Positions Offered`,
         `MD Senior Applicants`,
         `Total Applicants`,
         `MD Senior Matches`,
         `Total Matches`) %>% 
  pivot_longer(cols = -c(Specialty, Year, Class),
               names_to = 'Name',
               values_to = 'Value') %>% 
  arrange(Specialty, Class, Year)

long_table_absolute
long_table_percent <-  merged %>% 
  select(Specialty,
         Year,
         Class, 
         `MD Senior % Filled`,
         `Total % Filled`) %>% 
  pivot_longer(cols = -c(Specialty, Year, Class),
               names_to = 'Name',
               values_to = 'Value') %>% 
  arrange(Specialty, Class, Year)

long_table_percent 
p <- long_table_absolute %>% 
  filter(Specialty == 'Dermatology', Class == 'PGY2') %>% 
  ggplot(aes(x = Year,
             y = Value,
             color = Name)) +
  geom_point() +
  geom_line() + 
  facet_wrap(~Specialty + Class, scales = 'free_y') +
  theme_bw() +
  scale_color_carto_d(palette = 'Bold') 

p

input <- 'Dermatology'
plot_data <- merged %>% filter(Specialty == input,
                                Class == 'PGY2')#long_table_absolute %>% filter(Specialty == 'Dermatology')

fig <- plot_ly(plot_data, x = ~Year)
fig <- fig %>% add_trace(y = ~`Positions Offered`, name = 'Positions Offered', mode = 'lines+markers')
fig <- fig %>% add_trace(y = ~`MD Senior Applicants`, name = 'MD Senior Applicants', mode = 'lines+markers')
fig <- fig %>% add_trace(y = ~`Total Applicants`, name = 'Total Applicants', mode = 'lines+markers')

fig %>%
  layout(                        # all of layout's properties: /r/reference/#layout
    title = input, # layout's title: /r/reference/#layout-title
    xaxis = list(           # layout's xaxis is a named list. List of valid keys: /r/reference/#layout-xaxis
      title = "Time",      # xaxis's title: /r/reference/#layout-xaxis-title
      showgrid = F),       # xaxis's showgrid: /r/reference/#layout-xaxis-showgrid
    yaxis = list(           # layout's yaxis is a named list. List of valid keys: /r/reference/#layout-yaxis
      title = "uidx")     # yaxis's title: /r/reference/#layout-yaxis-title
  )
No trace type specified:
  Based on info supplied, a 'scatter' trace seems appropriate.
  Read more about this trace type -> https://plotly.com/r/reference/#scatter
No trace type specified:
  Based on info supplied, a 'scatter' trace seems appropriate.
  Read more about this trace type -> https://plotly.com/r/reference/#scatter
No trace type specified:
  Based on info supplied, a 'scatter' trace seems appropriate.
  Read more about this trace type -> https://plotly.com/r/reference/#scatter
No trace type specified:
  Based on info supplied, a 'scatter' trace seems appropriate.
  Read more about this trace type -> https://plotly.com/r/reference/#scatter
No trace type specified:
  Based on info supplied, a 'scatter' trace seems appropriate.
  Read more about this trace type -> https://plotly.com/r/reference/#scatter
No trace type specified:
  Based on info supplied, a 'scatter' trace seems appropriate.
  Read more about this trace type -> https://plotly.com/r/reference/#scatter
p <- long_table_absolute %>% 
  filter(Specialty == 'Dermatology') %>% 
  ggplot(aes(x = Year,
             y = Value,
             color = Name)) +
  geom_point() +
  geom_line() + 
  facet_wrap(~Specialty + Class, scales = 'free_y') +
  theme_bw() +
  scale_color_carto_d(palette = 'Bold') 

ggplotly(p)

Export

long_table_absolute %>% write_tsv('../../analysis/output/long_table_absolute.tsv')
long_table_percent %>% write_tsv('../../analysis/output/long_table_percent.tsv')
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKCmBgYHtyfQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShwbG90bHkpCmxpYnJhcnkocmNhcnRvY29sb3IpCmBgYAoKYGBge3J9Cm1lcmdlZCA8LSByZWFkX3RzdignLi4vb3V0cHV0L21lcmdlZF90YWJsZS50c3YnKQpgYGAKCmBgYHtyfQptZXJnZWQgJT4lIAogIGZpbHRlcihTcGVjaWFsdHkgPT0gJ0Rlcm1hdG9sb2d5JykgCmBgYAoKYGBge3J9CmxvbmdfdGFibGVfYWJzb2x1dGUgPC0gbWVyZ2VkICU+JSAKICBzZWxlY3QoU3BlY2lhbHR5LAogICAgICAgICBZZWFyLAogICAgICAgICBDbGFzcywKICAgICAgICAgYFBvc2l0aW9ucyBPZmZlcmVkYCwKICAgICAgICAgYE1EIFNlbmlvciBBcHBsaWNhbnRzYCwKICAgICAgICAgYFRvdGFsIEFwcGxpY2FudHNgLAogICAgICAgICBgTUQgU2VuaW9yIE1hdGNoZXNgLAogICAgICAgICBgVG90YWwgTWF0Y2hlc2ApICU+JSAKICBwaXZvdF9sb25nZXIoY29scyA9IC1jKFNwZWNpYWx0eSwgWWVhciwgQ2xhc3MpLAogICAgICAgICAgICAgICBuYW1lc190byA9ICdOYW1lJywKICAgICAgICAgICAgICAgdmFsdWVzX3RvID0gJ1ZhbHVlJykgJT4lIAogIGFycmFuZ2UoU3BlY2lhbHR5LCBDbGFzcywgWWVhcikKCmxvbmdfdGFibGVfYWJzb2x1dGUKYGBgCmBgYHtyfQpsb25nX3RhYmxlX3BlcmNlbnQgPC0gIG1lcmdlZCAlPiUgCiAgc2VsZWN0KFNwZWNpYWx0eSwKICAgICAgICAgWWVhciwKICAgICAgICAgQ2xhc3MsIAogICAgICAgICBgTUQgU2VuaW9yICUgRmlsbGVkYCwKICAgICAgICAgYFRvdGFsICUgRmlsbGVkYCkgJT4lIAogIHBpdm90X2xvbmdlcihjb2xzID0gLWMoU3BlY2lhbHR5LCBZZWFyLCBDbGFzcyksCiAgICAgICAgICAgICAgIG5hbWVzX3RvID0gJ05hbWUnLAogICAgICAgICAgICAgICB2YWx1ZXNfdG8gPSAnVmFsdWUnKSAlPiUgCiAgYXJyYW5nZShTcGVjaWFsdHksIENsYXNzLCBZZWFyKQoKbG9uZ190YWJsZV9wZXJjZW50IApgYGAKCmBgYHtyfQpwIDwtIGxvbmdfdGFibGVfYWJzb2x1dGUgJT4lIAogIGZpbHRlcihTcGVjaWFsdHkgPT0gJ0Rlcm1hdG9sb2d5JywgQ2xhc3MgPT0gJ1BHWTInKSAlPiUgCiAgZ2dwbG90KGFlcyh4ID0gWWVhciwKICAgICAgICAgICAgIHkgPSBWYWx1ZSwKICAgICAgICAgICAgIGNvbG9yID0gTmFtZSkpICsKICBnZW9tX3BvaW50KCkgKwogIGdlb21fbGluZSgpICsgCiAgZmFjZXRfd3JhcCh+U3BlY2lhbHR5ICsgQ2xhc3MsIHNjYWxlcyA9ICdmcmVlX3knKSArCiAgdGhlbWVfYncoKSArCiAgc2NhbGVfY29sb3JfY2FydG9fZChwYWxldHRlID0gJ0JvbGQnKSAKCnAKYGBgCgpgYGB7cn0KaW5wdXQgPC0gJ0Rlcm1hdG9sb2d5JwpwbG90X2RhdGEgPC0gbWVyZ2VkICU+JSBmaWx0ZXIoU3BlY2lhbHR5ID09IGlucHV0LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIENsYXNzID09ICdQR1kyJykjbG9uZ190YWJsZV9hYnNvbHV0ZSAlPiUgZmlsdGVyKFNwZWNpYWx0eSA9PSAnRGVybWF0b2xvZ3knKQoKZmlnIDwtIHBsb3RfbHkocGxvdF9kYXRhLCB4ID0gflllYXIpCmZpZyA8LSBmaWcgJT4lIGFkZF90cmFjZSh5ID0gfmBQb3NpdGlvbnMgT2ZmZXJlZGAsIG5hbWUgPSAnUG9zaXRpb25zIE9mZmVyZWQnLCBtb2RlID0gJ2xpbmVzK21hcmtlcnMnKQpmaWcgPC0gZmlnICU+JSBhZGRfdHJhY2UoeSA9IH5gTUQgU2VuaW9yIEFwcGxpY2FudHNgLCBuYW1lID0gJ01EIFNlbmlvciBBcHBsaWNhbnRzJywgbW9kZSA9ICdsaW5lcyttYXJrZXJzJykKZmlnIDwtIGZpZyAlPiUgYWRkX3RyYWNlKHkgPSB+YFRvdGFsIEFwcGxpY2FudHNgLCBuYW1lID0gJ1RvdGFsIEFwcGxpY2FudHMnLCBtb2RlID0gJ2xpbmVzK21hcmtlcnMnKQoKZmlnICU+JQogIGxheW91dCggICAgICAgICAgICAgICAgICAgICAgICAjIGFsbCBvZiBsYXlvdXQncyBwcm9wZXJ0aWVzOiAvci9yZWZlcmVuY2UvI2xheW91dAogICAgdGl0bGUgPSBpbnB1dCwgIyBsYXlvdXQncyB0aXRsZTogL3IvcmVmZXJlbmNlLyNsYXlvdXQtdGl0bGUKICAgIHhheGlzID0gbGlzdCggICAgICAgICAgICMgbGF5b3V0J3MgeGF4aXMgaXMgYSBuYW1lZCBsaXN0LiBMaXN0IG9mIHZhbGlkIGtleXM6IC9yL3JlZmVyZW5jZS8jbGF5b3V0LXhheGlzCiAgICAgIHRpdGxlID0gIlRpbWUiLCAgICAgICMgeGF4aXMncyB0aXRsZTogL3IvcmVmZXJlbmNlLyNsYXlvdXQteGF4aXMtdGl0bGUKICAgICAgc2hvd2dyaWQgPSBGKSwgICAgICAgIyB4YXhpcydzIHNob3dncmlkOiAvci9yZWZlcmVuY2UvI2xheW91dC14YXhpcy1zaG93Z3JpZAogICAgeWF4aXMgPSBsaXN0KCAgICAgICAgICAgIyBsYXlvdXQncyB5YXhpcyBpcyBhIG5hbWVkIGxpc3QuIExpc3Qgb2YgdmFsaWQga2V5czogL3IvcmVmZXJlbmNlLyNsYXlvdXQteWF4aXMKICAgICAgdGl0bGUgPSAidWlkeCIpICAgICAjIHlheGlzJ3MgdGl0bGU6IC9yL3JlZmVyZW5jZS8jbGF5b3V0LXlheGlzLXRpdGxlCiAgKQpgYGAKYGBge3J9CnAgPC0gbG9uZ190YWJsZV9hYnNvbHV0ZSAlPiUgCiAgZmlsdGVyKFNwZWNpYWx0eSA9PSAnRGVybWF0b2xvZ3knKSAlPiUgCiAgZ2dwbG90KGFlcyh4ID0gWWVhciwKICAgICAgICAgICAgIHkgPSBWYWx1ZSwKICAgICAgICAgICAgIGNvbG9yID0gTmFtZSkpICsKICBnZW9tX3BvaW50KCkgKwogIGdlb21fbGluZSgpICsgCiAgZmFjZXRfd3JhcCh+U3BlY2lhbHR5ICsgQ2xhc3MsIHNjYWxlcyA9ICdmcmVlX3knKSArCiAgdGhlbWVfYncoKSArCiAgc2NhbGVfY29sb3JfY2FydG9fZChwYWxldHRlID0gJ0JvbGQnKSAKCmdncGxvdGx5KHApCmBgYAoKIyMjIEV4cG9ydApgYGB7cn0KbG9uZ190YWJsZV9hYnNvbHV0ZSAlPiUgd3JpdGVfdHN2KCcuLi8uLi9hbmFseXNpcy9vdXRwdXQvbG9uZ190YWJsZV9hYnNvbHV0ZS50c3YnKQpsb25nX3RhYmxlX3BlcmNlbnQgJT4lIHdyaXRlX3RzdignLi4vLi4vYW5hbHlzaXMvb3V0cHV0L2xvbmdfdGFibGVfcGVyY2VudC50c3YnKQpgYGAKCg==